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Cn ' Abstract. In 1999, Griinbaum, Haine and Horozov defined a large family of 

commutative algebras of ordinary differential operators which have orthogonal 
3 , polynomials as eigenfunctions. These polynomials are mutually orthogonal 

with respect to a Laguerre-type weight distribution, thus providing solutions 

^\ ' to Krall's problem. In the present paper we give a new proof of their result 

which establishes a conjecture, concerning the explicit characterization of the 
dual commutative algebra of eigenvalues. In particular, for the Koornwinder's 
generalization of Laguerre polynomials, our approach yields an explicit set of 

■^^ ' generators for the whole algebra of differential operators. We also illustrate 

r \ ' how more general Sobolev-type orthogonal polynomials fit within this theory. 
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H ' 1. Introduction 

C^ , 

In 1929, Bochner [3 characterized the classical orthogonal polynomials (Jacobi, 

Laguerre, Hermite and Bessel) as the only ones that satisfy a differential equation 

of the form 

B (x,— \pn{x) = \nPn{x), U £ Nq, (1.1) 

where B{x,-^) — X]7=o ^ji^)l^ i^ ^ second-order differential operator with coeffi- 
cients bj (x) independent of the degree index n and with eigenvalue An independent 
of the variable x. In 1938, H. L. Krall [T7] posed the general problem to find all 
families of orthogonal polynomials {Pn(a;)}5^o which are eigenfunctions of a differ- 
ential operator B ~ X^ilo ^j(^) daj? of arbitrary order. In particular, he showed that 
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2 P. ILIEV 

the order m must be even and two years later [18] he completed the classification 
for operators of order 4, by adding three new families of orthogonal polynomials. 
During the last 70 years numerous families of orthogonal polynomials which are 
eigenfunctions of higher-order differential operators B were constructed. 

Inspired by the bispectral problem [5], Griinbaum and Haine discovered in [7] 
that the new families introduced by Krall can be obtained from the classical ones, 
by applying the Darboux transformation to the Jacobi recurrence operator. This 
brought new "bispectral" techniques in this old problem. Up to a linear change of 
the variable x, all known solutions of Krall's problem belong to one of the families 
below: 

[i] The Laguerre-type polynomials introduced by Griinbaum, Haine and Horo- 
zov [8] by applying the Darboux transformation to the recurrence operator 
for the Laguerre polynomials. These polynomials are orthogonal on [0, oo) 
with respect to a weight distribution of the form 

wix) = — -x"e"^ + y^ (-l)-'u,-(5^^'^(x), where a 6 Nq. 

r(a-t- 1) ^ ^-^ y J J \ J' 

finitely many jGNq 

[ii] The Jacobi-type polynomials introduced by Griinbaum and Yakimov [9] 
by applying the Darboux transformation to the recurrence operator for 
the Jacobi polynomials. These polynomials are orthogonal on [—1, 1] with 
respect to one of the following weight distributions: 

^'^ ' " 2"+^+ir(a + i)r(/? + 1) ^ ' ^ ^ ' 

+ H {-iyujS^^\x-l), whereaeNo, (3 > -1, 

finitely many JGNq 

or 



^,M = n^+13 + 2) _ ^ 



+ J2 {-iyu,s(^\x-i) 

finitely many jGNq 

+ Yl {-iyvjS(^^\x + l), wheiea,PeNo. 

finitely many j'^iNo 

Both papers use the general bispectral techniques developed in [2]. Moreover, 
the results in |8l [9] imply that to each of the families above, we can associate 
a "nontrivial" commutative algebra V of differential operators B such that ()1.1|) 
holds for some A„. By "nontrivial", we mean that the operators in the algebra 
are not obtained by taking constant coefficient polynomials of a fixed operator. 
The proofs in [8l [9] guarantee the existence of operators in the algebra V for every 
sufficiently large even order, thus raising the natural question of giving an explicit 
characterization of this algebra (and in particular, the operator of minimal order). 
Recently [TI] , we studied the algebra of differential operators for the Jacobi-type 
polynomials and used the results to obtain multivariate extensions of the family [ii] 
within the context of partial differential operators invariant under rotations. The 
main goal of the present paper is to explain how these techniques can be adapted 
for the Laguerre-type polynomials and to establish a conjecture in [8] concerning 
the explicit characterization of the isomorphic algebra A of eigenvalues A„. We also 
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illustrate how our approach leads to explicit formulas for the differential operators 
and how more general Sobolev-type orthogonal polynomials fit within the above 
theory. This leads to a uniform description of a large family of commutative algebras 
of ordinary differential operators which have Laguerre-type orthogonal polynomials 
as eigenfunctions, containing all known examples. 

A more detailed description of the results in the paper is as follows. Using the 
Darboux transformation, one can naturally associate to each of the Laguerre-type 
weights listed in [i] a certain Casorati determinant T{n) G C[n]. We denote by A 
the algebra of all A„ € C[n\ such that A„ — A„_i is divisible by T{n — 1) and we 
prove that A C A (i.e. for every Xn £ A there exists a differential operator B 
satisfying ()l.ip ). Next, we show that if the parameters Uj do not satisfy a certain 
algebraic relation, then A = A. In other words for generic parameters {mj}, all 
possible eigenvalues A„ in equation (jl.ip are in A. This fact was observed in several 
examples in [8 . A similar rule applies (generically) in the purely continuous version 
of the bispectral problem (with the discrete derivative An — A„_i replaced by the 
usual one), see 5^. However, our methods and arguments are very different from 
the ones used in the pioneering work of Duistermaat and Grtinbaum. Although 
we focus mainly on the Laguerre-type polynomials [i] , the results apply directly to 
more general Laguerre-type polynomials orthogonal with respect to a Sobolev-type 
inner product, which we discuss briefiy in the examples. 

The paper is organized as follows. In the next section we introduce the Laguerre 
polynomials and the basic notations which we use in rest of the paper. In Section 
|3] we briefly review the sequence of Darboux transformations from the recurrence 
operator of Laguerre polynomials which leads to the families listed in [i] and the 
construction of T{n). We introduce certain free parameters {/3j} which parametrize 
the Darboux transformations and which are in one to one correspondence with the 
parameters {uj} appearing in the weight distribution. In Section |4] we prove that 
A C A. In Section [5] we show that for generic parameters {/3j} (or equivalently, 
for generic {uj}) we have A — A. At the end of the paper we list three instructive 
examples. In the first example, we consider the Koornwinder's generalized Laguerre 
polynomials [14] which correspond to weights of the form 

w(x) = —- —x"e~^ + uoS(x), where a € Nq. 

i [a + I) 

The computation of the differential operator of minimal order has a long history 
and goes back to the works of H. L. Krall, A. M. Krall and L. L. Littlejohn [151 
UM [TTI [TBI [T^ . The explicit form of this operator for arbitrary a G N was obtained 
by J. Koekoek and R. Koekoek [12]. Our techniques applied to this case show that 
A ^ A and yield an explicit set of generators for the whole algebra V. In the second 
example, we consider in a detail the case when 

w(x) — — -x°'e^^ + UQ^ix) — uiJ'(x), where a G Nq. 

r(a + 1) 

We write an explicit formula for the algebraic equation which allows to determine 
whether uo and u\ are generic or not. In particular, when w(x) = e^^ + uq8{x) — 
uiS'{x) we show that yl is a proper subalgebra of A when uq and ui are not 
generic. In the third example, we explain how our results apply to Laguerre-type 
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polynomials orthogonal with respect to the Sobolev inner product 



(/(^),.9(a:)) 



1 



r(a + i) 



f{x)g{x)x^e--dx+[f{0)J'mA 



■ff(0)' 
5'(0) 



where a G Nq and A is a symmetric 2x2 matrix. 

2. Laguerre polynomials 



In this section we review briefly the basic properties of Laguerre polynomials. 
For n e No and a > — 1 the Laguerre polynomials {L"{x}}^^q are defined by 

r(a + n+l) 



m^) 



iFi{~n;a + l;x). 



(2.1) 



r(a + i)r(n + i) 

Here and in the rest of the paper, we use the standard gamma function r(a;) and 
iFi denotes the confluent hypergeometric function of the first kind 



iFi{a;c;x) 



^ (c)fcfc! ' 
fe=0 ^ ' 



\^ 



It will be convenient to set U^Ax) = for every j G N. The polynomials L'^{x) are 
mutually orthogonal on [0, oo) with respect to the gamma distribution x°'e~^dx. 
Explicitly, we have 

LZ{x)LZ,{x)x''e--dx = Ii^±^i±ll <?„,„. 

n! 

What is important for us is that the Laguerre polynomials are simultaneously eigen- 
functions of a second-order difference operator in the discrete variable n, which is 
independent of x, and of a differential operator in the continuous variable x, which 
is independent of n. More precisely, they solve the discrete-continuous bispectral 
problem 

-{n + 1)L:^+,{x) + {2n + a + l)Ll(x) - {n + a)Ll_^{x) = xi^(x) (2.2a) 

From the explicit formula (|2.ip one can easily deduce the following differential- 
difference relation 

±[L-{x) - Ll^,{x)] ^ ~Lf,_,{x). (2.3) 

By abuse of notation, we shall identify a function /„ defined on No with the 
corresponding semi-infinite vector [/o, /i, /2, ■ . ■]*. In particular, if we think of 
I/" (a:) as the semi- infinite vector 

[L^{x),LUx),mx),...f 

then the recurrence relations (j2.2ap can be written in the compact form 

J(")i^(x) = a;L^(x), (2.4) 

where J^°'' is the tridiagonal semi-infinite (Jacobi) matrix 

6o ao 

ci 6i ai 

C2 62 02 ' (2.5a) 



J(") = 
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with entries 

a„ = — (n + 1), bn = 2n + a + 1, c„ = — (n + a). (2.5b) 

3. DARBOUX TRANSFORMATIONS 

In this section, following [5], we describe the result of k successive Darboux 
transformations starting from J*-"' for a £ N and k < a. 

First we review briefly the lattice version of the Darboux transformation for 
bi-infinite matrices. Let j7o be a tridiagonal bi-infinite matrix with nonzero off- 
diagonal entries. If we factor j7o as a product of an upper-triangular and a lower- 
triangular matrices, then we can produce a new matrix by exchanging the factors. 
Following [3D], we shall refer to this operation as the Darboux transformation. If 
we iterate this process k times, we obtain a new tridiagonal bi-infinite matrix J^ as 
follows 

Jo = VoQo rv Ji = QoVo = PiQi rv • 

Jk-i = Qk-2Vk-2 = Vk-iQk-i r^J = Jk = Qk-iVk-i. (3.1) 

From (|3.ip it follows that 

JQ = QJo (3.2) 

and 

Jo* = VQ, (3.3) 

where V = Po'Pi ■ ■ ■ Vk-i and Q = Qk-iQk-2 ■ ■ ■ Qa- Similarly to the semi-infinite 
case, we shall identify a function /„ defined for n S Z with the corresponding bi- 
infinite vector [. . . , /_2, /-I, /o, /i, /2, ■ • • ]*■ Thus we can think of the bi-infinite 
matrices as linear transformations acting on functions /„ by matrix multiplication. 
Equations (|3.2p and (|3.3p imply that 

ker(Q) C ker(J(f) and Jo(ker(Q)) C ker(Q). (3.4) 



Conversely, one can show that if (|3.4p holds then there exists a tridiagonal matrix 
J which is obtained by a sequence of Darboux transformations as in (|3.ip . 

The lower-triangular matrix Q is uniquely determined by its kernel, up to a 
multiplication by a diagonal matrix on the left. More precisely, if {^n }jZo is a 
basis for ker(Q) and /„ is an arbitrary function then 

Qfn = .g„Wr„(V'i°\ Vi'\ . . . , i^i'-'lfn), (3.5) 

for an appropriate function 5„. We use Wr„ to denote the discrete Wronskian (or 
Casorati determinant): 

Wr„(5ii),g(2), . . . ,5^) = det(gW^.+Ji<,,,<,. 

Combining the above remarks, we see that the sequence of Darboux transformations 
(|3.1[) for Jo is characterized by choosing a basis for ker(Q) satisfying 

JoVi"^ = and JoV^i^'^ = V-rV"^^ for j = 1, . . . , fc - 1. (3.6) 

We want to apply the above construction with Jb = J^'"-* where J''-"' is an 
appropriate bi-infinite extension of the semi-infinite Jacobi matrix J^"' for Laguerre 
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polynomials. We define 



J 



(") - 



hU 



"0 



where 



-(n+e+1), 6^j = fe„+e = 2(n+e)+a+l, c^, 



-n+e 



-(n+e+a). 



Since a £ N we see that for e ^ Z, the off-diagonal entries a^ and c^ are nonzero. 
For j G {0, 1, . . . , fc — 1} we define 



/,iJ 



(-Ip /n + j 



(1 - a)j V J 



and (hl'^ 



(-Ip fn + a + j 
j\ \ a+j 



(3.7) 



It is easy to show that the functions {(/>'i+e} of the variable n G Z are linearly 
independent and satisfy 



J 



{a]n,0 



0, fori = 1,2 



J("V:;^, = (^:;^\ forj = i,...,fc-i, * = l,2. 

Thus, we can write the functions f/'n as a linear combination of (fif^^ as follows 



V'n ^ = ^il3j-l(t>nie + 7j-/0fj+J- 



;=o 



',(" 



If 7o = then ^h — (3o ^ 0, i.e. we can take -0; 



(0) 



1. One can check that J7i 



in p.ip coincides (up to a conjugation by a diagonal matrix) with JT"^" ^^. This 
means that the operator J = Jh can be obtained by a sequence of (fc — 1) Darboux 
transformations starting from J^""^^^ . Thus we can assume that 70 7^ 0, hence we 
can take 70 = 1. Since Q depends only on the space span{'0„ , ■(/)„ , . . . , ?/;„ }, 
and not on the choice of the specific basis, we can can take 7^ = for j > 0. Thus 
we shall consider a basis for ker(Q) of the form 



^i^'=Eft-i0iie 



1=0 






(3.8) 
Pk-i- We shall also normalize 



depending on a G N and k free parameters /3o, /3i, 
the matrix Q by taking g^ = 1 in (|3.5p . 

Now we consider the limit e — ;> 0. Since lime^o O-i+e = 0, it follows that the 
intertwining relation p.2p holds for the semi-infinite parts of the bi-infinite matrices 
J'o, J and Q. In other words, we have 



JQ = QJ 



(a) 



(3.9) 
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where J*-"' is the semi-infinite matrix in (|2.5|) . J is a similar semi- infinite Jacobi 
matrix 

bo flo 

Ci &i fli 

C2 &2 a2 



J 



and Q is a lower-triangular semi-infinite matrix, acting on functions /„ by 

Qfn = Wr„(V'i°), Vi'),. . . , Vi'=-'\/„), (3.10) 

where 



^i^^^Y.^^-^'^n'+^'n', 



(3.11) 



1=0 



with the convention that fj = when j < 0. 

One can easily write explicit formulas for the entries of the matrix J. Indeed, if 
we consider the k x k determinants defined by 

T(n)=Wr„(V^(°),^W,...,V'i'=-i)) (3.12) 

and 



p{n) 



^f 3 €ls 



ni-fc Vn-k 



^r^^ 
e-/^ 
e-/^ 



Vn-k 



(3.13) 



then using equations (|2.5p . (|3.2I) and (|3.9p one can deduce that 
. _ r(n-l) 



-(^ + 1), 



Pi'"^ + 1) / 1X P{^) 



r(7i) 
-f a + 

T(n) 



r(n- 1) 



-(n -I- a — fc). 



(3.14a) 
(3.14b) 
(3.14c) 



T{n - 1) ' 

In particular, from p.7p and the above formulas we see that the entries of the 
matrix J are well defined rational functions of n when 

T(n) 7^0 for n = -1,0,1,2,.... (3.15) 

We shall call the set of parameters /? = (/3o, /3i, . . . , Pk-i) admissible if ^^^^ holds. 
In the rest of the paper we shall work with admissible /3. One can show that 

ni=i(a-j) ^ 

which means that t(— 1) ^ is equivalent to ^o 7^ 0- 

Finally, we denote by L"''^(x) the polynomials defined by 

L^^^(x) =gL^(x) = Wr„(^(°),^W,...,Vi'=-i),L^(x)), (3.16) 

which depend on the free parameters a £ N, and j3 — (/3o, /3i, . . . , Pk-i)- 
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From ([231), dSS]), ([XTO]) and ([3l6| it follows that 

JL^-'^ix) = xL'^f^ix). (3.17) 

From equations p.l4ap and p.l4cp we see that the off-diagonal entries a„ and c„ 
of the matrix J are nonzero. Therefore, by Favard's theorem, there exists a unique 
(up to a multiplicative constant) moment functional A^ for which {L'^'^ {x)}'^^q is 
an orthogonal sequence, i.e. 

M{LZ'^{x)L'^'^{x)) = 0, for n j^ m and M ({Lf^'^{x))^) ^ 0. (3.18) 

More precisely, one can show that there exist constants uq, ui, . . . , Wfc_i such that 
the moment functional M. is given by the weight distribution 

fc-i 

^(^) ^ ^— rTy:^"-'e-- + ^(-l)^-«,<5«(x), (3.19) 

\a K). ,^^ 

where 5 is the Dirac delta function. The parameters {uj} correspond to a different 
parametrization of the sequence of Darboux transformations p.ip , see ^ Theorem 
2, p. 287]. We shall continue to work with the parameters {Pj} and we give the 
explicit connection between {/3j} and {uj} in the examples when fc = 1 and fc = 2. 



Remark 3.1. In the proofs of the main results (Theorem 14.11 and Theorem 15. 3p 
we shall only use the fact that the functions ^n' are polynomials of n (which 
follows immediately from equations p.7p and p. lip ), but not their explicit form. 
Thus all statements will hold if we use p.l6p and arbitrary polynomials ipn to 
define extensions of Laguerre polynomials. However, these polynomials will not 
(in general) satisfy the three-term recurrence relation p.l7p (hence they will not 
be orthogonal with respect to a moment functional). One can show that these 
polynomials will satisfy a higher-order recurrence relation. Indeed, note that if /„ 
is a polynomial in n of degree s, then [J'*-"'']*"^^/n-i-e = 0. Using this fact we see 
that if we define Q by p. 101) with polynomial functions ipn , then there exist s e N 
and a semi-infinite (2s -I- l)-band matrix J which is zero everywhere, except for s 
consecutive subdiagonals on either side of the main diagonal such that 

The above relation shows that the corresponding polynomials defined by (I3.16P will 
satisfy a higher-order recurrence relation analogous to (j3.17p . Interesting examples 
of such polynomials, satisfying a Sobolev-type orthogonality were studied by nu- 
merous authors, see for instance H] and the references therein. The higher-order 
recurrence relations can naturally be connected to theory of matrix- valued orthog- 
onal polynomials [6] and the Toda lattice [1] . We discuss the commutative algebra 
of differential operators for a specific family of Sobolev-type polynomials in the last 
section which extends the polynomials considered in |13j . 

4. The commutative algebras A and T> 

In this section we prove that the polynomials L'^'^{x) defined in equation p.l6p 
are eigenfunctions for all operators in a commutative algebra of differential opera- 
tors. 
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Let D denote the associative algebra of ordinary differential operators generated 
by the operators Di — -£; and D2 ~ x-^ — x-^, i.e. 

\ dx ' dx'^ dx I 

Thus S is an associative algebra generated by two elements satisfying the comniu- 
tativity relation {Dx.D-^ — D\ — D^. 

Note that equation (j2.2b[) can be rewritten as 

^ d 

BLZix) = nLl{x), where B = -xj^ - {a + 1 ~ x)— G S). (4.1) 

Recall that T{n) defined by (|3.12l) is a polynomial of n. We denote by A the 
algebra of all polynomials h[n) such that h{n) — h{n — 1) is divisible by r(n — 1) in 

C[n]: 

A = ihin) e CN : '^i^)-^i^-^) e CNJ . (4.2) 

I T(n-l) J 

It is not hard to see that A contains a polynomial of every degree greater than 
deg(r). We can define the elements of A as "discrete" integrals of polynomials 
divisible by T{n). To make this more precise we introduce the following notation: 
for m,n £ Jj and for a function /(s) defined on Z it will be convenient to use the 
following notation 

fELm+l/(5) ifn>TO 

/(s)d/id(s) ^ N if n = 771 

l-Er=„+i/(s) if'^<TO. 



Thus 



^n— 1 

/(s)(i/id(s) = f[n) + / f{s)d^d{s) for all m, n G Z. 



Note that ii h e A then 



and 



, , h(n) — h(n — I) ^, . ,.„, 

5('^) = t \, e C n (4.3) 

T[n — 1) 



pn pn—1 

h{n)^ g{s)T{s - l)dfid{s) + h{0) ^ g{s + l)T{s)dfid{s) + h{0). (4.4) 

With the above notations we can formulate the main result in this section. 

Theorem 4.1. For every h G A there exists Bh G S such that 

S,L^^^(.T) = h{n)L':-'^{x), (4.5) 

for all n G Nq. Thus, V ~ {B^ '■ h G ^1} is a commutative subalgebra ofT), 
isomorphic to A. 

For the proof of Theorem 14.11 we shall need two lemmas. First we formulate a 
discrete analog of a lemma due to Reach [21] . 
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Lemma 4.2. Let /„ , /„ , . . . , /A be functions of a discrete variable n. Fix 

ni, n2, . . . , Jik+i £ ^ o,nd let 

F^ = E(-l)'^'^'/«'^ /" /i°^Wr,(/(i), ... JO), ... , /('=+i))d^,(,), (4.6) 

wzi/i the usual convention that the terms with hats are omitted. Then 

Wr„(/(i), . . . , /«, K) = / /i°)Wr,(/ii), . . . , /i^-))dMd(.) 

Jn, + , (4.7) 

xWr„(,/W,..., /('=+!)). 

The above lemma was used in [ID i to construct commutative algebras of differ- 
ential operators for specific soliton solutions of the Toda lattice hierarchy. Since the 
proof is simple and we need it for one of the main results in the paper, we briefly 
sketch it below. 



Proof of Lemma \4-S\ Note that 



f(i) 



f(2) 



Jn Jn 

/•(I) f(2) 

Jn-1 Jn-1 

/•(I) /-(Z) 

■I n-k+l J n-k+l 
/•(I) f(2) 

Jn-l Jn-l 



Ak+l) 
Jn 

f(fe+l) 
Jn-l 



/•(fc+1) 
Jn-k+1 
Ak+l) 
Jn-l 



= 0, for every Z = 0, 1, . . . , fc — 1. 



Expanding the above determinant along the last row we obtain 

fc+i 

^(_l)fe+i+,_^W^Wr„(/(i), . ..f^\..., /('^+i)) = for ; = 0, . . . , fc - 1. (4.8) 

.3=1 

Using (|4.7p and (|4.8p we see that for Z = 0, . . . , fc — 1 we have 

k+l pri 

F„_,=^(-l)'=+i+Vi^, / /i°)Wr.(/W,..., /«,..., /f+i))dMd(.), (4.9a) 

„ -I J n^ 



and 



i=i 



fc+i 



F„_,=5:(-l)'=+i+Vi^J /i°)Wr.(/W,..., /«,..., /i'=+i))dMd(.) 

j=i -'". (4.9b) 

-/rWr„(/W,...,/r^)). 



If we plug (|4.9ap and (|4.9b|) in Wr„(/„ , . • ■ , /« ,Fn), then most of the terms cancel 



(fe) 



by column elimination and we obtain (|4.7I 



D 



Remark 4.3. We list below important corollaries from the proof of Lemma l4T2l 
(i) Note that the right-hand side of (|4.7p does not depend on the integers ni,n2 ■ ■ ■ ,nk- 
Moreover, if we change Uk+i then only the value of 



/i°)Wr,(/W,...,/W)d/id(.) 



"fc+i 
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will change by an additive constant, which is independent of n and /„ . Thus, 
instead of (|4.6I) we can write 

leaving the lower bounds of the integrals (sums) blank and we can fix them at 
the end appropriately. This would allow us to easily change the variable, without 
keeping track of the lower end. 
(ii) From (|4.8p it follows that for every I — — 1,0,1...,A; — 1 we can write F„ also as 

Fn = E(-l)'^'^Vi^'^ / /i"^Wr.(/ii), ...,/«,..., fi'+'^)d^,,is), 

and changing the variable in the discrete integral we obtain 

Fn = E(-1)'+'^V^-) / /i>r,(/« , . . . , /if,, . . . , /itrVA^d(s). (4.10) 

(iii) Let us apply (ii) with I — [^y^J and consider the sum consisting of the first k 
integrals: 

^' pn 

^« = ^(-D'^+i+.'/W / /(°) wr.(/W , . . . , /« , . . . , fiir')dMs). 

(k) 

Expanding each Wronskian determinant along the last column we can write Fn 
as a sum of k terms Fn '™ , each one involving as integrand one of the functions 
fs+m 1 where m = I ~k + l,l — k + 2, . . . ,1. We can use (|4.8I) once again, this time 
for the functions /„ ',..., fn (i.e. omitting /„ ), to change s as follows: 



• If TO > we can replace n with n — m in the upper limit of the integral, 
or equivalently, if we keep the upper limit of the integral to be n, we can 
replace s by s — to in the integrand. Thus Fn '"^ will have /i as 
integrand (and the integration goes up to n). 

• If TO < —1 we can replace s by s — m — 1, thus Fn will have JsLi as 
integrand (and the integration goes up to n) . 

If /n J /n 1 ■ ■ ■ , fn are polynomials of n, then summarizing the above observations 
we see that Fn can be written as a sum of terms of the form 

r^°Hn) / [r(i)(s)/i'^-+i)+r(2)(s)/i!t'^]dA^d(s), (4.11) 



where r^'^'{n), r^^'{n) and r'^^'(n) are some polynomials of n. 

Lemma 4.4. For every r(°)(n), r(^^(n), r'^^'(n) G C[n] there exists S G S such 
that 

rW(n) / [r^''>{s)L:{x)+r(^Hs)LU{^)W^{s)^BL'^{x). (4.12) 
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Proof. Using (|4T|) and ([231) "^^ see that 
\ ax J 



d 

dx 
d 
dx 



(L^(x)-L^_i(x))d/id(. 



Llix). 



Using (|4.ip once again we obtain (|4.12l) with 



B 



r^^\B) 1 



d 
dx 



„(2) 



(s+i; 



dx 



M 



{B) e 2). 



D 



Proof of Theorem \4-. 1\ Let h G A and let g{n) G C[n] be the polynomial defined in 
(j4.3p . We apply Lemma IT^ with 



We are now ready to present the proof of the main result in this section. 



f,^P^^i'-'^ioi-j^l,2,...,k 
fi'^+'^ ^ L'^ix). 

Using p.l2p . (I3.16P and (|4.4I) we see that with the above choice of functions /„ , 
the right-hand side of equation of (|4.7p is equal to {h{n) + c)L"'^(a;), where c is a 
constant independent of n and x. Thus, the proof of (j4.5p will follow from (I4.7P 
if we can show that, for an appropriate choice of the integers rij, there exists a 
differential operator B G 1) such that 



F,,^BL"^{x). 



(4.13) 



Indeed, if (|4.13p holds, then we can pull B in front of the Wronskian determinant 
on the left-hand side of (|4.7p which combined with p.l6p will establish (|4.5I) . We 
use now Remark 14.31 Note that the last term in the sum in equation (|4.10p is 
a polynomial of n times L'^{x) hence, using (|4.ip . this term can be rewritten as 
B"L'il{x) for some operator B" e 1). Thus it remains to prove that for the sum 
Fn consisting of the first k terms we have 

F^'^) = B"'LZ{x) for some B'" e 2). 

Since Fn is a sum of terms of the form given in (14. lip , we see that the proof 

D 



follows immediately from Lemma 14.41 

Remark 4.5. Using the explicit formulas for B in Lemma |4^ one can easily write 
exact formulas for the operator Bh in (|4.5p . Moreover these formulas are relatively 
simple for small values of k. For instance, when k = 1 we have 

Bh = g{B + 1) r A _ i") r{B) + h{B + 1), 
where g is defined by (|4.3p and B is the operator given in (j4.ip . 
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Remark 4.6. It is perhaps useful to note that the order of operator Bh constructed 
in Theorem 14.11 equals twice the degree of the polynomial h. One way to see this is 
as follows. From equations (j2.2p one can deduce that 



{n + a)L'^_,{x) = 


— X 


dx- (" + ^)dj 




hence for every j £ N we have 


{n + a^j + l),K_^{x] 


= 


i2 

'X-—-{0L^\ 


d 

dx 



LZix). 

This combined with definition p.l6p of L"'''(x) shows that there exist a polyno- 
mial 9{n) (which is independent of a;) and a differential operator _B in a; (which is 
independent of n), such that 

e{n)LZ'^{x)^BLl{x). (4.14) 

Multiplying (|4.5I) by 9{n) and using equations (14.11) and (I4.14p we deduce that 

{BhB - Bh{B))L'^^{x) = for aU n e No, 
hence B^B — Bh{B). This shows that 

order(i?,i) ^ order(/i(B)) = 2deg„(/i(n)). 

Remark 4.7. Note that equation (|2.ip can be used to define L'^{x) for all n e 
C \ (— N). Thus, we can think of L"{x) as a meromorphic function of n on C 
with simple poles at — N. For a G N, we can use p.l6p to obtain an extension of 
L"'°{x) which is meromorphic in C. By construction, these functions will satisfy 
the bi-infinite second-order difference equation 

JL::if,(x)=x4X(^). (4-15) 

where J7 is the bi-infinite matrix constructed in Section[2]and we interpret L'^'_!^^{x) 
as a bi-infinite vector. One can show that the function L'^'^{x) will also satisfy the 
differential equations (|4.5|) . Indeed, note that L^'^(x) is an entire function of x. 
Expanding L^'^{x) in a power series of x, we see that equation (|4.5p will hold if 
and only if the coefficients of x'^ on both sides are equal for every fc e Nq. Since 
these coefficients are polynomials of n, we see that (14. 5p holds for all n G C \ (— N) 
if and only if it holds for all n G Nq. Thus, by Theorem 14.11 equation (|4.5p holds 
for every n G C\ (— N). This combined with equation (|4.15p yields the bispectrality 
of the meromorphic (in n) function L'^'^{x). 

5. The commutative algebras A and f) 

Equations (|3.17p and (|4.5p establish the bispectral properties of the generalized 
Laguerre polynomials iv"'^(a;). More generally, we can consider the algebra A of 
all complex valued functions h(n) defined on No for which there exists a differential 
operator 

^'^ = E^'(-)^ (5-1) 

such that 

BhL^''>^{x) = h{n)L'^-^{x) for aU n G Nq. (5.2) 
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We denote by V the corresponding algebra of differential operators 

■D = {Bh: he A}. 

Griinbaum, Haine and Horozov proved that A contains polynomials h{n) for every 
sufficiently large degree and observed that A — A in several examples, see [H 
p. 295]. Theorem 14.11 shows that A C A, V C V and we prove below that A = A, 
V — V ior generic set of parameters (5 = {/3o, /3i, . . . , f3k-i)- 

First we list some simple properties of A and 2?. If we consider equation (15. 2p 
for n = 0,1,2,..., nih we see that bj{x) is a polynomial of x of degree at most 
j. Thus, V C C(^,a;^). If we compare the coefficients of x" on both sides of 
equation (j5.2p we see that h(n) is a polynomial in n of degree at most rn^. In 
particular, this means that we can think of ^ as a subalgebra of <C[n\. 

We denote by TZ{hi{n), h2{n)) the resultant of two polynomials hi(n) and h2{n). 
From (12. ip we see that 

which combined with (13.161) shows that 

^(0) = Wr„(^(°),^«,...,4'=-\(-l)^(;^+^)). (5.4) 

From equations ^^7^, ((XTT|) . ([XT^ and ^^ we see that T{n) and ^^^^^(0) are 
polynomials of n with coefficients which are polynomials of total degree at most 
k in the free parameters /3o,/3i, . . . ,/3fc-i- In particular, we can evaluate T(n) and 
^^" (0) for arbitrary n G C 

Definition 5.1. We say that an admissible set of parameters /3 — (/3o, /3i, . . . , (3k-i) 
is generic if 

in(/3):=7e(r(n),L^^^(0))^0. (5.5) 

Explicit formulas for 9\{(3) in the cases fc = 1 and fc = 2 are given in the next 
section (see equation (|6.ip for fc = 1 and (|6.4p for A; = 2). 

Remark 5.2. From the observations above it follows that $H(/3) is a polynomial 
in the parameters /3o,/3i, ■ • ■ ,Pk-i- It is not hard to see that this polynomial is 
not identically equal to 0. Indeed, let us assume that 5H(/3) is identically equal to 
zero. Then there must exist no G C which is a common root of T{n) and L"''^(0) 
for every admissible set of parameters /3. From ()3.7|) . (j3.11|) and ()5.4p it follows 
easily that L"'''(0) is independent of /3fe_i, hence no must also be independent of 
Pk-i- But then tiq rnust be a root of aZ which is independent of I3k~2- This 

shows that no must be a root of —q^ — — which is independent of Pk-d.- Continuing 

this way (alternating riji) and Z/^''^(0)) we see that no must be independent of 
(3o, /3i, . . . , Pk-i- But this is impossible since for every fixed no € C, T(no) is 
a polynomial of total degree k in the parameters /?o,/3i, • • . ,/3fc-i which is not 
identically equal to (for instance, the coefficient of /3q is the nonzero constant 

The main result in this section is the following theorem. 
Theorem 5.3. For generic (3 — (/3o, /3i, . . . , /3fc-i) we have A — A and V = 'D. 
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It is easy to see that when k = 1, every admissible /3o is generic and therefore 
A = A and V — T) (see Subsection 16. ip . However, one can show that in the very 
next simplest case (a = fc = 2) we have A'^ A and D ^T? when /3 is not generic, 
see Subsection 16.21 

For the proof of the above theorem we need two auxiliary facts. 

Lemma 5.4. Let m e N and let f^^^n), f'-'^^n), . . . , /'•"'■^'^^n) be functions of n. 
Then 

Wr„(/(i)(n), /(2)(n), . . . , /('"+^)(")) Wr„(/(2)(n _ 1), /(3)(„ - 1), ... , /(™+i)(n - 1)) 

= Wr„(/(2)(n - 1), /(3)(n - 1), ... , /("+2)(n - 1)) Wr„(/(i)(n), f^^Hn), . . . , f^^+'\n)) 

- Wr„(/(2)(„), /(3)(„), . . . , /(™+2)(„)) Wr„(/(i)(n - 1), /^^H" " 1), ■ • ■ , f'^"'+'\n - 1)). 

(5.6) 

Proof. The proof of (j5.6p follows immediately by applying the Desnanot-Jacobi 
identity to the determinant Wr„(/(i) (n), /(2) (n), . . . , /(™+2) („)). Q 

Lemma 5.5. Let m G N and /et /'^-'^^(n), /(^^(n), . . . , /^™+-'^)(n) be polynomials of 
n. For j = m anrf j = m + 1, set G'-J\n) = Wr„(/(i)(n), /(2)(n), . . . , /(j)(n)). 
T/ien 

7e(G(")(n),G(™+i)(n)) ^ TZ (^G'^"'\n - 1) , G^"'+^\n)^ . (5.7) 

Proof. We prove the statement by induction on m. When tti = 1 it is easy to see 
that 

7^(G(l)(n),G(2)(r^)) = 7^ (G^^^n - 1), G(2)(n)) 

= n (/(i)(n - 1), /(i)(n)) 7e (/(i)(n), /(2)(„)) . 

Suppose now that the statement is true for some m E N and we want to prove it 
for m + 1 . For j = m + 1 and j = to + 2 let us denote 

gO)(n) = Wr„(/(2)(n), /(3)(n), . . . , /0)(n)). 
By the induction hypothesis we have 

7e(G(™+i)(n),g('"+i)(n)) = TZ (^G'^"'+^\n) , Q'^"'+^\n - l)y (5.8) 

From Lemma 15.41 it follows that 

■Jl ('g'("+i)("), G("+2)(n)) 7e (g^^+i^J^), Q^'^+^H'^ - 1)) 
= 7e (g("+i) (n), -Q(™+2) (n)G(™+i) (n - 1)) 

= 7e (G("+i)(n), g("+2)(n)) 7e (G("+i)(n - 1), G(™+i)(n)) , (5.9a) 

and 

n (G^"'+^\n - l),G(™+2)(n)) 7e ('G("+i)(n - l),g('"+i)(n - 1)) 

^ n (g(™+i)(" - 1), Q^^+^H" - l)G("+i)(n)) 

= 7e(G("+i)(n),g("+')(n)) 7e (G(™+i)(n - 1), G("+i)(n)) . (5.9b) 

The proof now follows from equations (|5.8I) and (|5.9p . D 
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Proof of Theorem \5.3\ We need to show that for generic [3, A C A. Let h be an 
arbitrary element of A. By definition, there exists a differential operator B^ G V 
such that (|5.2p holds. We know that Bh has polynomial (in x) coefficients and 
h{n) S <C[n\. Evaluating ()5.2|) at a; = and using the resulting equation for n and 
n — 1 we find 



LZ^fMiBnLZ'H^)) „ - K^^mBHL^^^ix)) 



x=n 



(5.10) 



= L:L^i(0)L^^^(0)(Mn)-Mn-l)) 

We look at both sides of ()5.10p as polynomials of n and we would like to prove 
next that r(ri — 1) divides the left-hand side. For fixed j G N we apply Lemma [5T4l 
with m = A: and fW(n) = ^(0), /(2)(n) = 7/;i"\ f(^^n) = t/;,^, . . . , /('=+i)(n) = 
^(^fe-i)^ /('=+2)(„) ^ i^(o) and we obtain: 

Wr„ (^Vf , V'i'\ . . . , V'r'\ ^(0),L^(0)j r(n - 1) 

= i^^r (0)^(0) -L-^(0)^(0). 

The last equation shows that T{n~l) divides L"i5i(0)^^^(0) -£°'^(0) '''^5'' (0) 
for every j S N which proves that T(n — 1) divides the left-hand side of (|5.10p . 
Therefore, r(n — 1) must divide the right-hand side of (|5.10p . It remains to show 
now that for generic /3 the polynomials T{n — 1) and L^_j^(0)I/"'^(0) are relatively 
prime. From Lemma 15.51 it follows that 

n (r(n), L^^^(O)) = n (rin - 1), L^^'^(O)) . 

Thus using the notations in Definition 15.11 we find 

7^ (r(n - 1), 4^L^i(0)L^^''(0)) = MP)f + for generic /3, 

completing the proof. D 

6. Explicit examples 
6.1. One Darboux step. Suppose that a G N and fc = 1. Then 



Tin) = V4°^ = /3o + 



n + a 



and formula p.l6p gives 

L^^(^«)(x)=Wr„(^(°),i^(. 

The polynomials LZ (x) wee orthogonal on [0, oo) with respect to the weight 
distribution 

w{x) = r^x'^^^e^^ + uoS(x), where uq — -—. 

(a-1)! /^o 

From the explicit formula above it is easy to see that 

^n -f a — 1 
a-1 



L^'(^°)(0) = -/3o(' 
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$H(/3o) =7^(r(n),L-(''o)(o)) ^ (-1)" J°^" ' (6.1) 



Combining the last formula with the formula for r(n) and definition Definition l5.ll 
we see that 

n2c 

\a Po 

Recall that /3o 7^ since we work with an admissible set of parameters (i.e. T{n — 
1) ^ for n G No) and therefore every admissible /3o is generic. This means that 
when fc = 1 we always have A = A and V = f). We write below an explicit set of 
generators for both algebras. 

For j G No we consider the polynomials of n defined by 

.»,„,. ,„(";-^).('';^^)(:-;j). (.., 

Note that 

r(n - 1) V ] J 

This shows that /i(-')(n) G ^forallj G No- Moreover, since deg„(/i^-'-'(n)) = a+j+1 
we see that 

A = <C[h^''\n),h^^\n),...,h^'^\n)]. 

Using Remark 14.51 we can write explicit formulas for the operators Bj := Bf^u): 

Bj = g^^\B + 1)(^- l\ t{B) + h^^\B + 1), 

where B is the second-order differential operator for Laguerre polynomials given 
in ()4.ip . The algebra of differential operators T) is therefore generated by the 
operators -Bo, Bi, . . . , Ba of orders 2a + 2, 2a + 4, . . . , 4a + 2. The construction of 
the operator Bq of minimal order 2a + 2 has a long history. For a = 1, a = 2 
and a = 3 the differential equation of minimal order goes back to the works of 
H. L. Krall, A. M. Krah and L. L. Littlejohn [HI [ISl [HI [H [19] . The explicit form 
of Bq for arbitrary a G N was obtained by J. Koekoek and R. Koekoek |12) . 

6.2. Two Darboux steps. Let us take now fc = 2 and let a be an integer greater 
than 1. The polynomials L°''^(x) = LT/'^"''^'\x) are defined by ((XTCl) with 

^r=/?o0r + 0?.^°=/3o+(";!^") (6.3a) 

^(1) = /3i<^r + Ml^' + ^Y = /3i + ^(n + 1) - f " + " + ^) . (6.3b) 

a— 1 \a+l/ 

They are orthogonal on [0, 00) with respect to the weight distribution 

(a -2)! 

where 

(a-l)(/?o + /3i) , a-1 

uq = -2 and ui = — - — . 

/?o /3o 

For L" (0) we obtain the following formula: 

L-(^-^^'(0) = -^(n + l)„_.[(" + ;^-l)+/3o 
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Using the explicit formulas for T{n) and L„ (0) one can derive the following 

formula for the resultant ?l(/3): 



OQ o5q — 5 



^a I 2 2/3oa J ^_ ^« + l 



a 



{a - 1)2"-2(q, + l)"[(a - 1)!]4" 

(6.4) 
Ignoring the unessential constant factors and /?q" ^ (which is nonzero) , we see that 
admissible /3o and /3i are generic if and only if they do not satisfy the algebraic 
equation 

Q-l 

(2a)"a!/3o"+i + [] [{a" -a- 2aj)f3o - (a' - l)/3i] = 0. (6.5) 

j=o 

Using the explicit connection between (mo,wi) and (/3o,/3i), one can easily rewrite 
the above equation in terms of uq and ui. 

Note that T{n) is a polynomial in n of degree 2a. Thus, the commutative algebra 
A is generated by polynomials of degrees 2a + 1, 2a + 2, . . . , 4q: + 1. The corre- 
sponding algebra of differential operators V is therefore generated by operators of 
orders 4a + 2, 4q; + 4, . . . , 8a + 2. 

Let us consider in details the simplest case a = 2. The weight distribution is 

w{x) — e~^ + uo6{x) — uiS'{x). 
Equation (|6.5p becomes 

or equivalently, in terms of uq = — (/?o + Pi)/ Po f^nd ui = l//3o we get 

9uq + ISuqUi + 5u^ + 32ui = 0. 

If /3o and /3i do not satisfy (|6.6p , then ^ = ^ is an algebra generated by polynomials 
of degrees 5, 6, . . . , 9 and P = I? is an algebra of differential operators generated 
by operators of orders 10, 11, . . . , 18. The operator of minimal order 10 is given 
explicitly in [51 pages 292-295] (in slightly different notations). 

However, if /3o and /3i satisfy (|6.6p . then one can show that A contains the 
following polynomial of degree 4 

h{n) =n^ + 271^ + (28/^0 - l)^^ + (36/3i + 28/?o - 2)n. 

Thus if (|6.6p holds we have A'^ A^T) '^V and the polynomials Un ' {x) satisfy 

a differential equation of order 8. 

6.3. Sobolev-type orthogonality. Let a be an integer greater than 1 and let 
us consider a two-dimensional subspace of the space span{(/)^'°, ip'^f' , 0^'^, 0^'"'^} with 
basis {V'n 7 '0n }• If we assume that the coefficient of ^^-^ in ■0n and the coefficient 
of (^2.1 in i/^n are nonzero, we can pick ■0„ and ij^n as follows: 

V^i") = /3o0i'° + ;o0^-^ + 0^'° = /?o + -^(n + 1) + f" + ") (6.7a) 

a — 1 \ a J 

^W=A0i^V;i0^'^-f0^-i=/3i + -^(n + l)-f" + ';+^). (6.7b) 

a— 1 \a+l/ 
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If we take lo — and h = /3o equations (|6.7|) reduce to (|6.3p . If we define a bi- 
infinite matrix Q with basis {V'„+e, V'n+e} then using the notations in Section [3] we 
have 

ker(Q) C ker[j(")]2, 

but for arbitrary Iq and h 

j(")(ker(Q))^ker(Q). 

Since only the first part of equation (|3.4I) holds, we cannot expect to find a tridiago- 
nal matrix J', which is obtained by a sequence of two Darboux transformations from 
J7'"', but we can find a pentadiagonal matrix J which is a Darboux transformation 
of [j(")]2 as follows: 

[j.(a)]2 =-pQrxJ =QV. 

In the semi- infinite limit (e — )> 0), this construction leads to a semi- infinite upper- 
triangular matrix P and a semi-infinite lower-triangular matrix Q, with two diag- 
onals above (resp. below) the main diagonal: 



P 



Po,o 



Poa Po,2 

Pl,l Pi, 2 Pi, 3 
P2,2 P2,3 



P2,4 



90,0 

9i,o 9i,i 

92,0 92,1 



92,2 



satisfying 



[ j(")]2 ^ PQ and J = QP, 



(6.8) 



where J is a semi-infinite pentadiagonal matrix. 

Let us consider the polynomials L" "' (a;) — Ln{x) obtained by applying 

Q to L"{x), or equivalently by using the Wronskian formula p.l6p with k = 2 and 
V"" , V'n given in (I5T71) . i.e. 

L^{x) = QL-{x) = YfT^{^k^^\^^\Ll{x)). (6.9) 

As we noted in Remark 13. 1[ Theorem 14.11 and Theorem 15.31 can be applied for 
I/„(x) (since the proof used only the fact that ipn' are polynomials of n). Thus, 
if we define r(n) and A by equations (|3.12p and (|4.2p . then for every h ^ A there 
exists Bh G £>, such that 

BhLn{x) = h{n)Ln{x). 

We can also consider the algebras A, T) defined at the beginning of Section [5] and 
for generic parameters /Sq, l3i,lo, li we will have A = A and V = V. 

The polynomials Ln{x) will not (in general) be orthogonal with respect to a 
moment functional. However, if we pick "carefully" the free parameters /3o, /3i , ^o, h 
the polynomials Ln{x) will be orthogonal with respect to the Sobolev inner product 



{F{x),G{x)) = 



1 



(«-2)!yo 
where A is a symmetric 2x2 matrix 

A = 



F{x)G{x)x°'-^e-''dx + [F{0), F'(0)] A 






G(0)" 
G'(0)J 
(6.10a) 



(6.10b) 
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Indeed, note that equation ()6.9p gives 

Ln{x) = Qn.nL'^ix) + (7„,„_iL"_i (x) + q„,„_2L"_2 (a:) , 

which combined with the orthogonahty relations for L"{x) shows that for the inner 
product defined in (|6.10p we have 

{Ln{x), x^) —Q when 2 < j < n. 

Thus, we need to pick the parameters so that (Z„(x),x^) = when j G {0, 1} and 
n > j. From equations (I6.8p . (I6.9P and (|2.4p we see that PLn{x) — x^L'^{x), or 
equivalently 

Pn,nLn{x) + _p„,„+lI/„+l (x) + Pn.,n+2Ln+2{x) = x"^ L"(x) . 

It is easy to see that the numbers p„^,i-f 2 are nonzero, and therefore the last equation 
shows by induction that 

• if {Ln{x), 1) = for n =: 1 and n = 2, then {Ln{x), 1) = for all n e N. 

• if {Ln{x),x) — ior n — 2 and n — 3, then {Ln{x), x) = for all n >2. 

Summarizing the above observations we see that the polynomials Ln{x) will be mu- 
tually orthogonal with respect to the inner product (|6.10l) if and only if (L„(x), 1) — 
for n G {1, 2} and (L„(x), x) = for n G {2, 3}. If det(A) 7^ 0, one can show that 
these four equations are satisfied if we pick f]o,PiJo,li as follows 

{a~l)ui (a - l){uo + ui) ^rii.n 

^" = -^^aT' ^' = — d^m — ' ^ ^ ^ 

{a ~ l)«o {a - l){ui + vq) /«11K^ 

^° = ^^aT' ^' ^ d^^) — ■ ^^-''^^ 

If Vq = 0, then Iq ~ 0, li = j3o and we obtain the polynomials in the previous 
example (iteration of two Darboux transformations). For vq ^ 0, the polynomials 
are no longer associated with a moment functional (since (x,x) ^ (x^, 1)). 

Since deg„(T(n)) = 2a the commutative algebra A is generated by polynomials 
of degrees 2a + 1, 2a + 2, . . . , 4a + 1. In the case when the matrix A is diagonal 
(i.e. wi = 0) the operator of minimal order 4a + 2 was computed explicitly in |13) . 

Finally, if the matrix A is singular, we can still apply the above construction by 
a limiting procedure, but we need to pick a different basis for the two-dimensional 
subspace of span{0^'",0^''',0^'-'^,(/)^'-'^}. For instance, if we want to consider the 
inner product (|6.10ap with 

"0 0" 
vo\ ' 

then the polynomials _£„ (x) will be given by (|6.9p with ipn and ipn defined by 

V'i°^=n+1 (6.12a) 

^W=(" + :|V(" + ")-^^. (6.12b) 

Va+ly \ a J vo 

Note that in this case deg„(T(n)) — a+1, which means that our theory will produce 
a differential operator in V of minimal order 2a -I- 4, in agreement with the results 
in [13]. 
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